F(R) gravity equation of state 
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1. Introduction. — In a key article 1'], Jacobson ob- 
tained Einstein's equations (Ee) from local thermo- 
dynamics arguments only. To derive this result, he 
generalized black hole (Bh) thermodynamics to space- 
time thermodynamics as seen by a local observer. He 
noted his finding strongly suggests that, in a fun- 
damental context, Ee are to be viewed as an equa- 
tion of state and, therefore, they should probably 
not be taken as basic for quantizing gravity. This is 
consistent with the idea that gravity is an emergent 
phenomenon of a more fundamental framework, like 
string theory (e.g. [4]). Were this true, not only Gen- 
eral Relativity, but presumably all generalized gravity 
theories should be seen under this same light. 

Modified gravity models constitute a very impor- 
tant dynamical alternative to ACDM cosmology, in 
that they have the capability to describe the current 
accelerated expansion of our Universe (dark energy 
epoch), but also the initial de Sitter phase and infla- 
tion, and even the galaxy rotation curves correspond- 
ing to dark (and ordinary) matter [3|. We will here 
prove that Jacobson's derivations can be generalized 
to cover these more complicated theories of gravity 
that are extensively used nowadays. First, we review 
Jacobson's arguments to introduce the basic notions 
and then derive the desired generalization. In par- 
ticular, we will completely close the program for the 
so called f(i?) gravities (see e.g. [3|, and references 
therein), where the Lagrangian only depends on the 
Ricci scalar and its covariant derivatives, leaving the 
problem open for more general cases. In [J], the field 
eq. for f (i?) of polynomial form were derived using 
non-equilibrium thermodynamics arguments (see also 
[5|). Here, we propose an alternative approach where 
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local thermodynamic equilibrium is maintained, using 
the idea of "local-boost-invariance" introduced in [61]. 
2. Jacobson's construction in brief. — Any free-falling 
local observer p has some gauge freedom to describe 
his local coordinate system. The equivalence principle 
can be used to describe space-time in a vicinity of p 
as flat. Then, we choose the local space-like area ele- 
ment perpendicular to the world- line of p to have zero 
expansion rate 9 and shear cr at a given point on the 
history of p, that we call po- In this setting, the past 
horizon of po is called the "local Rindler horizon" at 
Po- Since, locally, we have Poincare symmetry, there 
is an approximate Killing field K generating boost at 
Poi vanishing at po, which we take future pointing to 
the inside past of pq. 

Having this basic setting, we are ready to give pre- 
cise meaning to the local thermodynamic definitions. 
First, note that local Rindler horizons are null and 
act as causal barriers. Therefore, we can associate 
entropy S* to it, measuring the "many degrees of free- 
dom outside" , what presumably results in entangle- 
ment entropy just at the horizon. With this under- 
standing, entropy is proportional to the area elements 
of the horizon, where a fundamental length has to be 
provided to give an uv cut-off. Heat Q is energy flow 
of microscopic degrees of freedom across the causal 
barrier, and is felt, therefore, via gravitational energy, 
where its source is undetectable. Lastly, the local tem- 
perature T is defined as "Unruh temperature" , as seen 
by a local accelerated observer hovering just inside the 
horizon. Energy flow has to be measured by this same 
observer, for consistency. 

In more detail, different accelerated observers would 
measure different energy flows and temperature, both 
diverging at the horizon but with constant ratio, and 
this is just what will be used. We have also imposed 
^ = a = at P07 to give a sort of "local definition of 
equilibrium" since, in general, causal horizons change 
in time as they expand and twist. In this construe- 



tion, locally and at the instant po there is not such a 
deformation, and the space is "at equilibrium" . 

2.1. Accelerated observer and approximations. — On 
the above point po with its associated local Rindler 
horizon 7i, take an accelerated observer hovering just 
inside the horizon x- By the above construction, x 
is an approximate local boost Killing field future di- 
rected to the past of po- Then, the variation of heat 
(caused by energy flow across the horizon), measured 
by Xj is SQ — J^TabX^'^"': where the integration is 
over a pencil of generators of H a,t pq. If iiT is a 
tangent vector to the generators of H, with affine 
parameter A such that A = at po, we have that 
X° = —k\K°' + O(A^), where k is the acceleration of 
X. Therefore, dE" = —k\ dX dA, where dA is the cross 
section area element of H. Thus, the final expression 
for the variation of heat, at leading order, is 



SQ = - kXTabK'^K'dXdA. 
In 



(1) 



Note that for x the Unruh temperature T is set to be 



T = k/2Ti . 



(2) 



On the other hand, Jacobson uses that the variation of 
the entropy is proportional to the variation of the hori- 
zon area A^ i.e. 5S — rjSA^ with 77 an unknown pro- 
portionality constant. Here, 5A measures the change 
of the area as we approach the point poj and therefore 
is given by 5 A = J^9{X)dXdA. Next, we use Ray- 
chaudhuri's equation to integrate 9 near po- In this 
coordinate system, at leading order in A, we obtain 
= -XRabK^K^ + 0{X^), and the relevant expres- 
sion for the entropy variation to this order is 



5S 



ri I XRabK''K''dXdA. 



(3) 



n 



2.2. Thermodynamic relations. — To derive informa- 
tion of thermodynamic systems, like the equation of 
state, we need just the basic thermodynamic relation 



5Q = T5S, 



(4) 



identities for the Ricci tensor, we get the integrability 
conditions: -^rjWaR = —'^ ah, and therefore the final 
form of the thermodynamic relation is 



Mah — nab — -Tr9ab 

V J \ 2 



Affab, 



(6) 



where A is an integration constant. 

To summarize, we have here obtained the Ee as an 
equation of state for a local free-falling observer. To 
deduce the above, we have used the following critical 
assumptions: (i) Measurements are done in a vicinity 
of a general point pq. (ii) Our local coordinate system 
is at equilibrium, in the sense that 6, a — at po. 
(iii) The accelerated observer x tends to i^, a null 
vector generator of the causal horizon, (iv) We always 
restrict ourselves to the leading order approximation 
in the afline parameter A. 

3. The general case of modified gravity. — We apply the 
above construction to more general theories of grav- 
ity. Following Iyer and Wald [6|, we just assume that 
our Lagrangian is diffeomorphism invariant, in a n- 
dimensional oriented manifold M, being the dynam- 
ical fields a Lorentz signature metric gab and other 
matter fields ip. The most general Lagrangian is 

L = ^\gab,Rcdef,^aiRcdef,...,'^(ai---'^a„)Rcdef, 
V',Va,l/',...,V(,, ...V,„)V'). (7) 

The corresponding field equations con be found by a 
variational procedure on {gab, V')i so that we get 



^L = e (E/'^ffafc + Ev,(5V) + d@, 



(8) 



where e is the volume element and a (n — 1)- 
form. Hence, the field equations of the theory are: 
E °^ — 0, E^ — 0. In 6], it was found how to write 
them from a variation of the {gab, Rcdef), as if they 
were independent variables, so that we get, after the 
corresponding identifications. 



xn ab 

^9 



< 



rp pqra p b 



2'^pVqEj^ 



pabq 



(9) 



and the functional dependence of S with respect to 
the energy and size of the system. In our case we have 
Eqs. (|1I2I3P at disposal, to get the beautiful relation 



T^bK^K" 



27r 



T^RabK^KK 



(5) 



Since K is an arbitrary null vector on Ti, we can write 
the unprojected equation Tab — -^rjRab + gabh, with 
h an unknown function, arbitrary as of now. Using 
then that the Ihs is divergence- free, plus the Bianchi 



where {Ag'', Ej^"" ') are the variations of L with re- 
spect to {gab,Rpabq) in each case, taken as indepen- 
dent variables. In the above expressions, if the deriva- 
tives of Rcdef occur in the Lagrangian, one integrates 
by parts and then takes its variation, to obtain E^"^ ''. 
This form of the field equations is useful due to 
its relation to Bh thermodynamics. Basically, it has 
been known for a while now 7] , that in the case when 
we have a stationary Bh solution, the entropy S can 
be calculated as a Noether charge evaluated at the 



bifurcation (n — 2)-surface of the event-horizon E. In 
these cases, the entropy is given by 



5* 



-27r / E 



1 ahpq 



ab^pqi 



(10) 



where eat is the binormal vector of E. 

What is less understood is the case of dynamical Bh 
entropy. There, the event-horizon is not bifurcated 
and, therefore, the above formula does not hold. Nev- 
ertheless, in [6|] a prescription that passes the basic 
tests of consistency for the corresponding entropy is 
presented, although it is not completely clear if that 
is a good answer. In any case, the idea is to approx- 
imate the metric 5, in a vicinity of a given point p 
of the event-horizon, by a boost-invariant metric g^"^. 
This is done by altering the original Taylor expansion 
of the metric around p, so that the new metric is boost 
invariant up to some order g, that defines the size of 
the vicinity where our approximation is valid. Then, 
for this boost invariant metric, there is a Killing vec- 
tor field that, on the horizon, is null, and vanishes at 
p. We thus have created an approximated bifurcation 
surface of order q and can use the same expression 
as before for the entropy, only that the integration is 
done on the boost-invariant variables: 



^dyn(^p) 



-2tt 



E 



ahpq 



^ab^pqi 



(u) 



ahpq 
R 



^i? (5^') (see [g for details). 



where E 

Having understood these modifications for calculat- 
ing the Bh entropy, we are almost ready to continue. 
Still, some information on the geometrical meaning of 
Eq. (fTO|l . In 8] it was noticed that, for a static Bh, 
entropy can always be re-expressed as the area of the 
bifurcation {n — 2)-surface A divided by 4 in units of 
an effective Newton constant Geff, i.e. 



S = 



A 



4G 



where 



1 



eff 



SttG 



eff 



T-i abpq 



(12) 



The above result has been checked for some string the- 
ory cases where it was found that Geff is indeed con- 
stant on the bifurcation surface. This result has to be 
supplemented with the key observation that the above 
effective Newton constant plays also the role of an ef- 
fective gravitational coupling for graviton exchange. 
In other words, the kinetic term of the n-dimensional 
graviton, obtained from the general Lagrangian L, is 
precisely of the form \E^ ^' {yrhbq^^hap + . ■ •) and, 
hence, E^ ^' can be thought of as the strength of 
the graviton interaction in all possible polarizations. 
In retrospective, Geff corresponds to the strength of 



the gravitational interaction in the particular polar- 
izations defined by the binormal of the bifurcation 
(n — 2)-surface A. 

3.1. Field equation as equation of state. — Now that 
we have the basic inputs for the possible geometri- 
cal interpretation of the Bh entropy S for generalized 
theories of gravity, we are ready to consider the prob- 
lem of defining there a local version of Bh thermody- 
namics, following the steps of Jacobson. Note that 
all definitions regarding local observers (accelerated 
or not), local Rindler horizons, and so on, are based 
on differential geometry and the equivalence principle. 
We expect all these definitions to hold in the general 
setting and, thus, we leave them unchanged. 

What makes the difference, being one of the key 
points in this generalization, is the definition of the 
local entropy. After the discussion above, it seems 
natural to relate entropy to the area of the causal 
horizon, only that now we replace the proportionality 
constant with a field-dependent effective constant. In 
other words, we state that the local variation of the 
entropy is still proportional to the variation of the 
area of the causal horizon, but in units of this effective 
Newton constant. Therefore, we write now [10| 



6S = 6iT^eA), 



(13) 



where rje is, in general, a function of the metric and 
its derivatives to a given order 1 + 2, i.e. 



Ve = Ve 



igab, Rcdef , ^ Rp 



(14) 



Using the above ansatz, we are ready to proceed with 
our derivation. Since we just change the definition of 
entropy variation, due to the energy flow across the 
local Rindler horizon, we get the modified expression 



5S 



X {l]e Rab - VaVbTleWK'^dX dA + 0{\^ 



H 



(15) 

It is important to notice that, in this expression, 
{rie,k°''^ aVe) src to be evaluated at its leading con- 
tribution in A. We have used its boost-invariant part 
at first order in lambda to effectively incorporate the 
boost invariant notion of [6| creating an "approxi- 
mated bifurcation point at first order in A" at po- 

The other part of the derivation is unaffected and 
gives the same result of ([T]), namely 



5Q 



( k\TabK''K^d\dA + 0{\^). (16) 

JH 



Therefore, the thermodynamic relation ([4]) implies 



TabK-'K" = l^ileRab ~ V ^a^ b)Ve)K^ K" . (17) 



At this point we consider the general differential equa- 
tion, removing the contraction with K, thus 



T 



ab _ Ve_ 
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j^ab _y{ayb)V^^^ab^^ ^^g) 



27r 



where the new terms are added based on the fact that 
K isa. tangent vector of the null geodetics at poj gener- 
ating local boost. Hence, at this point we have a local 
equation with two unknown functions of the metric 
and its derivatives {r]e,H). 

To find the form of these three functions, we use the 
integrability condition 

Va (^R'''' - V^'^ V^) ^ + g'^'^H) = 0, (19) 

V27r zTT / 

obtained from the observation that the rhs should be 
divergence free. After some algebra. 
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(20) 



At this point, with no lose of generality, we can set 

(21) 



"-"^^T.- 



finally obtaining the expression 



'7e V76i 



= 2:-\7''R + V''h. 

4tt 



(22) 



3.2. The specific case of f(R) gravities. — Eq. p9)) can 
in principle be solved in many different ways. Here we 
will consider the simplest possibility, that eventually 
leads to the so-called f(i?) gravities [J], a special — 
but phenomenologically very important — family of 
the general theories of modified gravity where only 
the Ricci scalar is involved, as well as their covariant 
derivatives, i.e. L = f (i?, V"i?). 

The solution of Eq. (|19p we are considering is per- 
haps the simplest one where the first two terms cancel 
each other. The last one can be easily integrated as- 
suming r/e is a function of R only, thus h — —-^{{R), 
with 2^ — ■§j^- Therefore, the final form once we have 
collected all terms in the above derivation is 



p.ab ^ ^f^ab _ y{a^b)^ 

dR OR 



V 



^2 5f 

V — 
dR 



(23) 
f a function of R and its covariant derivatives only. 

Eq. ([23]) is in fact the correct field theory equation 
for f (i?) gravities, provided we identify the function 
f(R) as the Lagrangian of the theory. Also, in this 



case the effective Newton constant of p2)) is related 
to T^e, as is expected from the relation 

i — pPV-s - _('„P'' gs _ gr psN 



dt ^ Ve 

OR " 2tt 



(24) 



Note that, for these theories, the different polariza- 
tions of the gravitons only enter in the definition of 
the effective Newton constant through the metric it- 
self. This is an important simplification that, in turn, 
permits to find the solution of the integrability con- 
dition (|19p . In retrospective, and to summarize, we 
have succeeded in our thermodynamic derivation of 
f(-R) gravities where, remarkably, exactly as in the 
case of Einstein gravity, the local field equations can 
be thought of as an equation of state. 

It will be very interesting to see if this derivation can 
be extended to the more complicated cases, stemming 
from string theory, where the full Riemann tensor is 
involved in the Lagrangian. This seems to imply a 
sort of tetrad decomposition of the effective Newton 
constant such that one recovers, at the end, only the 
polarization normal to the causal barrier of the local 
Rindler horizon. Work along this line is in progress. 
As a last comment, in our derivation we have used the 
first law, but no information is given bout the second 
law. In fact it is not know if the second law is present 
in generalized gravities (see [9])- 
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